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CD . ABSTRACT 

This is an exploratory study of the hghtest non-singlet scalar qq state on the lattice with 
^ ■ two dynamical quarks. Domain wall fermions are used for both sea and valence quarks on a 
^ . 16^ X 32 lattice with an inverse lattice spacing of 1.7 GeV. We extract the scalar meson mass 

1.58 ± 0.34 GeV from the exponential time- dependence of the dynamical correlators with 
^vai = ^sea and Nf = 2. Since this statistical error-bar from dynamical correlators is rather 
large, we analyze also the partially quenched lattice correlators with m^ai T^^sea- They are 
positive for m^ai > f^sea and negative for m^ai < f^sea- In order to understand this striking 
effect of partial quenching, we derive the scalar correlator within the Partially Quenched 
ChPT and find it describes lattice correlators well. The leading unphysical contribution in 
Partially Quenched ChPT comes from the exchange of the two pseudoscalar fields and is also 
positive for m^ai ^^sea and negative for m^^/ <?^sea at large t. After the subtraction of this 
unphysical contribution from the partially quenched lattice correlators, the correlators are 
positive and exponentially falling. The resulting scalar meson mass 1.51 ±0.19 GeV from the 
partially quenched correlators is consistent with the dynamical result and has appreciably 
smaller error-bar. 



1 Introduction 



The interest in the hght scalar mesons has been renewed recently [1] . The existence of the 
scalar mesons above 1 GeV is well estabhshed experimentally and there are enough scalar 
states between 1 GeV and 2 GeV to represent the scalar qq nonet [2]. The excess of one 
observed state in this region has been suggested as an indication for the glucball [3]. The 
lightest iso-triplct state above 1 GeV is ao(1450). The only scalar states below 1 GeV, which 
are experimentally well established, are iso-triplet ao(980) and iso-singlet /o(980) [2]. The 
existence of the complete scalar qq nonet roughly below 1 GeV would require another iso- 
singlet and two strange iso-doublets. The experimental evidence for the existence of a broad 
iso-singlet a meson around 600 MeV is growing [2,4], while the existence of the strange 
iso-doublet k, reported in [5] is even more controversial at present. This raises a question 
whether the lightest scalar qq states lie below 1 GeV or above 1 GeV. In the latter case, the 
observed scalar states below 1 GeV have to be interpreted as exotic states like qqqq [6], tttt 
or KK molecules, etc. 

In this paper we address the determination of the mass of the lightest scalar qq state 
with non-singlet flavor structure (referred to as the meson [2]), the long-term goal being 
to relate this state to the observed resonance ao(1450) or ao(980). We determine the mass 
of the gq meson using a lattice simulation of dynamical QCD (m^a/ = '^'sea) and partially 
quenched QCD {rrivai 7^ iTisea) with Nf — 2 degenerate sea quarks in both cases. Since our 
aim is qq state composed of the light u and d quarks, we employ Domain Wall Fermion 
(DWF) formalism, which has good chiral properties [7]. We comment also on the mass of 
the su and sd scalar mesons. 

While DWF formulation ought to be helpful in the long run, at present our numerical 
work has serious limitations. We are working with two dynamical flavors, which is not full 
QCD. Furthermore we have results only at one lattice spacing on a lattice box that is not 
very large and also quarks are relatively heavy. For these reasons this is an exploratory work. 
These issues can of course be improved with more computing resources. 

Before we introduce our work, we briefly review the recent lattice simulations of the 
light non-singlet scalar states. We quote only the statistical error-bars on masses since the 
continuum and infinite-volume extrapolations have not been performed in these simulations: 

• Fully quenched simulations of qq: 

The quenched qq correlator in the chiral limit was simulated by Bardeen et al. [8, 
9] with Wilson fermions. The correlators were found to be negative at small quark 
masses, which was attributed to the similar mechanism as observed in the present 
partially quenched study. The effects of quenching were modelled using the Quenched 
Chiral Perturbation Theory and subtracted in order to extract the scalar meson mass 
m„o = 1.326(86) GeV [9]. 

The RBC Collaboration simulated non-singlet and singlet scalar qq states with Domain 
Wall Fermions [10]. The quenching effect, which leads to negative correlators at small 
quark masses, was subtracted as in [8,9]. The result is rUaQ — 1.43(10) GeV if only 
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the leading chiral loop (one bubble) is taken into account, and m.ao — 1.04(7) GeV if 
next-to leading chiral corrections are included by resummation^. 

Mixing of the glueball and qq was studied in [3]. The quark mass was around and 
no attempt was made to go to the chiral limit. 

• Dynamical simulations of qq: 

The SCALAR Collaboration made an extensive simulation of the singlet qq state [11] 
and extracted also the mass of the non-singlct state to be niao ~ 1.8 GeV at rriT^/mp ~ 
0.7. They consider this estimate as an upper bound on the mass since they fitted the 
correlators at relatively low times, where contribution of the excited state might be 
sizable. 

UKQCD extracted niao ~ 1-0(2) GeV from the dynamical and the partially quenched 

simulation of qq [12]. Since they simulated only rriyai > f^sea, they did not observe the 
striking effect of partial quenching discussed below. For this reason they were able to 
extract the scalar mass from the exponential timc-dcpcndcncc. 

MILC [13] simulated qq state with three dynamical flavors and saw an indication for 
the intermediate state nrj, since this state is lighter than oq state at the lightest quark 
masses. 

• Alford and Jaffe reported an indication for the bound singlet and octet qqqq states 
below 1 GeV [6]. 

All simulations above employed Wilson fermions, except for RBC and MILC simulations, 
which employed Domain Wall and staggered fermions, respectively. 

The only simulation which employed chiral fermions to study light scalar mesons is the 
quenched simulation of RBC [10]. Chiral symmetry is expected to be particularly important 
for the singlet scalar meson a, which is intimately connected with the chiral symmetry 
breaking. Good understanding of the non-singlet correlator in the chiral limit is the first 
step toward the controlled study of the cr meson. As already mentioned, the present paper 
presents the dynamical simulation (m^ea = "^^a/) of the non-singlet qq correlator with Domain 
Wall fermions. We also simulate partially quenched QCD with rriyai 7^ ^sea- Two degenerate 
sea quarks have the range of masses corresponding to 500 — 700 MeV [14]. The scalar 
correlators for rusea — 0.02 and rusea — 0.03 at various rriyai are shown in Fig. 1. The 
correlators for rriyai ^ ''^sea are positive and have more or less exponential timc-dcpcndcnce. 
On the other hand, the correlators for rUyai < rUsea are negative due to a striking effect of 
partial quenching. We note that the point-point correlator should be positive definite in the 
dynamical QCD based on unitarity, which is broken in the partially quenched QCD. We 
derive the effect of partial quenching on the scalar correlator using the Partially Quenched 
Chiral Perturbation Theory (PQChPT) [15] in a finite box. The leading unphysical effect is 
due to the exchange of the two pscudoscalar fields, it is represented by the bubble diagram 

^The higher order chiral corrections C'(M^/(47r/)-^) have smaller effect on the scalar mass in Ref. [8,9] 
since their value of the pseudoscalar decay constant / is larger than the physical value. 
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in Fig. 3b and has no unknown parameters. We show that the bubble diagram gives 
a positive contribution for m^ai > "^sea and a negative contribution for m^ai < ^sea at 
large time separations. We find that the negative lattice correlators with m^ai < ^sea are 
well described by the bubble contribution. This enables us to extract mao in the partially 
quenched simulation. 
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The remainder of this paper is organized as follows. The details about the lattice simu- 
lation are presented in section 2. The dynamical correlators are analyzed in section 3. The 
resulting error on the scalar mass is rather large, which motivates us to analyze also the 
partially quenched correlators. The partially quenched artifacts on the scalar correlator are 
derived within PQChPT in section 4 and used to analyze the partially quenched correlators 
in section 5. Section 6 summarizes the conclusions on the mass of Oq meson and briefly 
comments on mass of the k resonance, while section 7 summarizes the general conclusions. 



2 Numerical simulation 

The RBC Collaboration has undertaken a large-scale simulation with Nf = 2 flavors of 
dynamical Domain Wall quarks with degenerate masses [14]. This is an improvement over 
the quenched simulations and represents an important step toward the simulation of QCD 
with three dynamical quarks of physical masses. The scalar correlators were calculated on 
the dynamical configurations with the volume N^Nt = 16^32 and a single lattice spacing, 
so we will not be able to extrapolate the scalar mass to the continuum and to the infinite 
volume in the present work. The configurations were generated using DBW2 gauge action 
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[16]^ with (3 = 0.80 and Domain Wall fermion action [7] with Mg = 1.8 and = 12 
{Lg is the extent in 5th dimension) [14]. The separate evolutions were performed for three 
different bare sea-quark masses m^ea = 0.02, 0.03, 0.04, which correspond approximately to 
M^~500— 700 MeV. The measurements of the correlators were performed on configurations 
separated by 50 HMC trajectories. Dynamical Domain Wall fermions have good chiral 
properties even at finite Lg with the additive shift in the mass due to the residual chiral 
symmetry breaking, rrires, being approximately 0.0014 [14], much smaller than either the 
input sea or valence quark masses. The inverse lattice spacing was determined from the 
p-meson mass and the preliminary result is ~ 1.7 GeV [14]. The current uncertainty of 
the lattice spacing has a small effect on the scalar mass; the uncertainty of the scalar meson 
mass in the present work is dominated by the statistical errors of the scalar correlators. 

The scalar correlators were measured for the degenerate valence quark masses mi = m2 = 
rriyai in the range rriyai = 0.01 — 0.05. These valence quark masses correspond approximately 
to M;r~380 — 770 MeV (Fig. 4). We simulated the correlators with p = 0, point source and 
point sink via 

^ J](0|g(f,t)rg(a;',t) q{y,0)Tq{y,0)\0) 

^ x,x',y 

with r = 7 on the lattices with un-fixed gauge. The average over the configurations with 
un-fixed gauge gives the point-point correlator [19] 

= ]^E^O|^(^'^)r^(^'^) g(y,0)rg(y,0)|0) . (1) 

This method enabled us to calculate also the singlet scalar and singlet pseudoscalar correla- 
tors and to determine the hairpin insertion itlq. 

The summary of the scalar and pseudoscalar correlators analyzed in the present paper is 
given in the Table 1. 

3 Analysis of dynamical correlators with m^ai = ^sea 

The mass iriaQ and the unrenormalized decay constant fao can be extracted from the dy- 
namical scalar correlators using the exponential fit in the conventional way. Indeed, we 
will verify that the additional contribution from the exchange of two pseudoscalar fields in 
PQChPT (Fig. 3b) exactly vanishes for m^ai — T^sea-, = 2 and mo — > oo (12), so the 
simple exponential fit is well justified. The extracted masses and decay constants are shown 

^DBW2 gauge action is known to break the reflection positivity of the transfer matrix [17], which is the 

counterpart of the unitarity for Euclidean lattice theory. Hence, there is a possibility that the negative value 
of the scalar correlator is partially caused by the complex eigenvalues of the transfer matrix. We concentrate 
our analysis at large time separation t > 6, where the negative contribution to the two-point function is 
expected to disappear. We perform the comparison between the full QCD and the (partially) quenched QCD 
with niyai < fnsea using the same DBW2 gauge action and find negative scalar correlator only in case of the 
(partially) quenched QCD. 
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Figure 2: The asterisks represent the mass niao ^-nd decay constant fao in the lattice units at the 
dynamical point m^ai = t^sea, which are obtained from the exponential fit of the scalar correlators. 
The dashed line and the value at ruq = are obtained with the linear fit. 



in Fig. 2, while Table 2 presents also the time ranges t = tmin — tmax and of the fit'^. 
The uncorrelated fits are used throughout this work and the error-bars are obtained using 
the jack-knife method. The linear extrapolation to the chiral limit m^ai = ^sea —>■ gives 
results in lattice units^ 

ruao = 0.93(20) , fao = 0.049(20) , (2) 

where the jack-knife error-bars are calculated as described in Appendix B of [18]. 

The resulting errors are rather large, which motivates us to extract the mass also from the 
partially quenched data with m^ai 7^ ^sea- This forces us to understand the effect of partial 
quenching in the following sections. The use of the Partially Quenched ChPT is crucial for 
this purpose since it enables us to subtract the significant partially quenched artifacts from 
the negative lattice correlators in case of rrivai < ^sea- 



4 Scalar correlator in partially quenched ChPT 

In this section we derive the non-singlet scalar correlator in the Partially Quenched ChPT 
(PQChPT) within the so-called p-expansion regime. We consider the theory with Njjai 
valence quarks (which can have different masses nii) and Nf degenerate sea quarks qs 
of mass rrisea- The theory incorporates also N^ai valence ghost-quarks of mass rrii, which 
cancel the closed valence-quark loops. PQChPT enables us to study of the partially quenched 
artifacts, which arise if the valence and the sea quark masses are not equal and if Nf ^ N^ai- 

^tmin is taken throughout this work high enough so that there is no visible effect of the excited states 
and at the same time as low as possible in order to avoid large statistical errors. The choice of tmax has 
negligible effect on the result and we take it at the time slice just before the signal is lost. 

"^The difference between the chiral extrapolations rUq and niq —rrires is negligible due to the 
smallness of rrires ~ 0.0014 [14]. 



6 



Our few lowest quark masses are low enough that M^/(47r/)^ <^ 1, while they are still large 
enough that M^L ^ 1 and we do not enter e-regime on our lattice. 




(a) (b) 

Figure 3: The contributions to the non-singlet scalar correlator in PQChPT: (a) The exchange 

of the scalar meson ao; (b) The bubble diagram is responsible for the unphysical effect of partial 
quenching and represents the exchange of two pseudoscalar fields The intermediate pseu- 

doscalar fields and can have the fiavor structure ~ qiqj, qiqj, (jiqj, ijiqs, qsQi , where 

qij are the valence quarks and qs the sea quark. 

Non-physical contributions to the scalar correlator in PQChPT arise from the exchange 
of pseudoscalar fields between qq source and qq sink. The leading contribution in the chiral 
expansion comes from the exchange of two pseudoscalar fields and is represented by the 
so-called bubble diagram in Fig. 3b. The two pseudoscalar fields $ and can be mesons 
^ij ~ QiQj and ~ (jiqs with Boson statistics, or mesons ~ q^qj with fermionic 
statistics. We do not consider the next-to-leading chiral corrections in PQChPT, which are 
suppressed by C(M^/(47r/)^) in comparison to the bubble diagram^. The lattice correlators 
can be interpreted as a sum of the ao-exchange at the tree level in Fig. 3a and the bubble 
diagram in Fig. 3b. 

For our purpose, we need the strong interactions of pseudoscalar fields in PQChPT [15] 
as well as the kinetic and the mass term for the ao field ^ 

£ = ^ stT[d^'Ud^U^ + /Vo stT[M^U + U^M] - ^(str + d^'aod^al - ml^aoal (3) 

with the physical values / ~ 95 MeV and three-flavor mo ~ 600 MeV-1000 MeV [20]. 
The field U = exp[\/2i$//] incorporates the SU{N^ai + Nf\N^ai)L x SU{N^ai + Nf\N^ai)R 
Goldstone field matrix $. The quark mass matrix A4 is 

M = diag{mi, mN,,i, rnsea, ■ ■ ■,msea , mi, ... , mjv^^J 

and the supertrace is defined as strA = Xla^a^a with = 1 for quarks and = —1 for 
ghost quarks. The parameter represents the slope of versus rUg and the pseudoscalar 
masses M are given by 

= 2^o(ma + mfe) , a,b^i,i,S, i = l,..,N^ai, S = l,..,Nf. (4) 

'""The NLO chiral corrections were taken into account for fully quenched scalar correlator by resummation 
[8 10]. They are more complicated in partially quenched theory because partially quenched theory implies 
several intermediate states in Fig. 3b, while quenched theory implies a single intermediate state irr]' . 

^Similar Lagrangian was used for the case of fully quenched ChPT in [8-10]. There /xq is denoted by ^vq. 
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As is standard, we neglect the Q;(9str term in the Lagrangian since a seems to be small 
[9,20]. We also need the couphng of Oq field and the pseudoscalar fields to the non-singlet 
scalar current^ 

mi ~ -/"o/'(f/ + f/^)i2 - \fvm /io /ao ao (5) 
where /ao plays the role of the scalar meson decay constant with a given normalization: 

(0|^2?i|ao) = - Vl28 //o /ao ■ 
The point-point scalar correlator with external momentum p = (0, E) 

C^^'it) = J2{^\q2{x,t)q,(x,t) ?-i(0,0)g2(0,0)|0) (6) 

X 

is computed on the lattice and can be related to the prediction of PQChPT 



aO 



I continuum 



+ B{t) . (7) 



The first term is the conventional Oo-exchange, while the second term B{t) = F.T.[B{p)] 
represents the contribution of the bubble diagram in Fig. 3b. The lattice Euclidean momenta 
p — 2sin(pE/2) and the discrete Fourier Transform F.T. on are used when we compare 
PQChPT predictions (7) with the lattice correlators^. The bubble diagram is calculated 
from the Lagrangian (3) and the current (5) in the Appendix A, giving 




(8) 



1 1 

k'^ + Mls r 1 1 2 

N~f {k + pf + Ml^ 1 , k'+^is [(FTm^ ^ (FTm|)2 + (fc2 + M2J(A;2 + M|2) 

Here k denotes the momenta in the loop (Fig. 3b) and the sum is performed over the 
allowed loop-momenta in the finite box of the lattice. We use the lattice Euchdean momenta 

k = 2sin(A;^/2) and k + p = 2sin((A;£' +pe)/2) when comparing PQChPT prediction (7) 
with the lattice correlators. The detailed expression for C^^'-"^^^{t) used in case of our finite 
lattice is presented in Appendix B. 

We note that the contribution of the bubble diagram (8) has no unknown parameters, 
since one can fix the values of M, /iq and mo from other considerations. We determine 

''Similar current was used for the case of fully quenched ChPT in [8]. There /Uo is denoted by ^ro- 
^When we refer to conventional "exponential fit" , we extract niao and fao from the fit to the first term 
in (7) using p = 2sin(j3^;/2) and the discrete F.T. 
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the pseudoscalar meson masses M and fXo — M^/(4mg) from pion correlators on the same 
lattices. The hairpin insertion mo~600 MeV— 1000 McV (normalized for 3-fiavors) has been 
determined from the rj' correlator in a number of references [20], but the exact value of itlq 
is not essential for the present work since the extracted scalar mass is almost independent 
of mo in the wide range mo = [600 MeV, oo] as will be demonstrated below. 

In order to understand the effect of the bubble contribution on the lattice correlators, we 
derive the asymptotic form of B{t) at large t for a correlator with p = on a lattice with 
a — > 0, qNt oo and finite qNl. The asymptotic form for degenerate valence quarks with 



mass rrivai is 



TD(j.\ t^°9 2/io 



' + n...2 ^. + (9) 



M'vv 3M2^ I (M2 - Ml,yf + - M^^ 

where M^, = Mgg + ^Nfml denotes rj' mass in a theory with Nf flavors. We have assumed 
Mr,' > Mvv in derivation of (9), which is satisfied for the pseudoscalar masses of physical 
interest. The asymptotic behavior is dominated by a pair of zero-momentum pseudoscalar 
fields with mass Myy = 4yUom.ya/ , a pair with mass Mvs = '^^^o{f^vai+^sea)i oi a pair with 
masses Myy and Af^/. The dominant contribution at large t for m^ai > f^sea is proportional 
to e~'^^^^^ and has positive sign. The dominant contribution for m^ai < ^sea is proportional 
to i-e"^^*^^* and has negative sign given by Myy — Mj^. The bubble contribution is inversely 
proportional to the spatial volume of the lattice, so the effect of the bubble contribution is 
much less important for larger lattices. We summarize these findings which apply^ for any 
values of Nf and mo as follows: 

The scalar correlator receives a positive contribution e^™""* from the exchange of 
meson and an additional bubble contribution from the exchange of two pseudoscalar fields 
$i$2- The bubble contribution is proportional to e^^^i+^^^'/A^I at large t and it is positive 
for m^ai ^ "^sea ^-nd negative for m^ai < ^sea- The scalar correlator with rriyai < rUsea 
has a negative sign at large t if the bubble contribution dominates over the ao-^xchange. 
The bubble contribution is particularly important on lattices with smaller spatial volume 
if Ml + M2 < m^Q. The magnitude of the bubble contribution is predicted by PQChPT 
(8), which enables us to extract mao and fao by fitting the scalar correlators to PQChPT 
prediction (7). 

We close this section by demonstrating that the analytical expressions for the bubble 
contribution (8,9) reduce to the known expressions in the fully quenched and the fully un- 
quenched limits: 

• The fully quenched ChPT corresponds to the limit rusea — > 00 or equivalently to the 
hmit Nf — > 0. Eqs. (8) and (9) reduce in this hmit to 



^These findings apply as long as the condition M^' > Myv is satisfied, which was assumed in the 
derivation of (9). 
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(p) = -16,? E (^^^ . (10) 

9, ,2 2 p-2Mt 

with M = Mvv for degenerate valence quarks. Expression (10) agrees with the fully 
quenched expressions used in [8-10].^° 

In the case oi full ChPT with SU{Nf) flavor symmetry mi — m.2 — msea, the expressions 
(8) and (9) reduce to 



,,2 



■ -{M+M ,)t -2Mt 

4 ^ + (7V2 4)^— 



for general uiq , (13) 



where M = Myv — Mss- This is in agreement with the conventional ChPT result. 
Bose symmetry and conservation of isospin allow only one intermediate state nr]' in 
Nf = 2 ChPT and this intermediate state gives vanishing contribution in the mo — > oo 
limit. There are three intermediate states nrj, KK and ttt;' in Nf — 3 ChPT. 



5 Analysis of lattice correlators using PQChPT 

Our scalar correlators (6) are calculated in the partially quenched lattice QCD for Nf — 2 
degenerate sea-quarks with mass m^ea — 0.02 — 0.04 and degenerate valence quarks with 
mass rriyai = 0.01 — 0.05 (Table 1). In this section we fit the lattice correlators to the one- 
loop prediction of the PQChPT given in Eq. (7). The PQChPT prediction is the sum of 
the ao-exchange diagram and the bubble diagram (Fig. 3). The magnitude of the bubble 
contribution (8) is completely determined for a given choice of mo since we determine the 
values of M and Hq from the pseudoscalar correlators on the same configurations: 

• The pseudoscalar masses are given by M^^ = |(M^^ + M^^), where Maa is determined 
from the pseudoscalar correlator with rriyai = rua and the sea quark mass of interest. 
The from our lattice correlators are listed in Table 3 and shown in Fig. 4. 

• We fix /io = M^-ai vai / {'^''^^vai) foT glvcu m,,,ai and rriseaj where nires is neglected since 
it is much smaller than either rrivai or nisea- Here rrivai is the input bare mass of the 
valence quark. Myai,vai is the mass of the pion with two valence quarks of mass rrivai 
at nisea of interest. 

^°The mg/S can be viewed as the hairpin insertion in the quenched theory with one valence flavor, while 
Ho = ^ro in [8-10]. 
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Figure 4: M'^irriyai) at rusea = 0-02 in lattice units. 



The bubble diagram vanishes in the case of the dynamical theory {rUvai 



nisea) with 

Nf — 2 and TOq — (12). In this case, the lattice correlator is interpreted solely by 
the ao-exchange (7), it has exponential time-dependence and the corresponding niao was 
extracted in Section 3. The bubble diagram is non-zero in general, so it has to be taken into 
account when the lattice correlators arc fitted by equation (7) in order to extract rriao and 
fao- The bubble contribution incorporates the physical contributions from the exchange of 
two pseudoscalars and also the unphysical effects of partial quenching when niyai 7^ ^sea- 
We note that our scalar correlators and the quark masses, used to determine //q = 
/ {Anig) , are not renormalized. We would like to emphasize that this does not prevent us 
from extracting the physical mass niao- This can be seen by rearranging Eqs. (6) and (7): 



Y,{0\q2qi{x,t) qiq2{0,0)\0) = coef. 



-niaoiNT-t) 



+ F.T 



MtB{p) 



16 



• (14) 



The product of quark mass and the scalar current m^gg is invariant under renormalization. 
The second term on the RHS of (14) depends only on the hadron masses and mo, which 
are also invariant under renormalization. The scalar mass mao can be therefore extracted 
from the first term on RHS without ambiguity from the renormalization. 



5.1 Analysis of dynamical correlators taking into account -kt]' in- 
termediate state 

The two-flavor dynamical scalar correlator receives a contribution from the exchange of oq 
and from the exchange of ttt]' . The contribution from vrr/' state vanishes in the limit tjiq oo 
and the fit reduces to the standard exponential fit used in section 3. The contribution of vrr^' 
intermediate state at finite itlq is given by the bubble contribution (8) with m^ai = f^sea and 
Nf = 2. The fit of the dynamical correlators to PQChPT prediction (7) at various mo gives 
mao and fao in Fig. 5. The results are almost independent of mo and agree with the result 
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of the conventional exponential fit in section 3. For this reason we will refer to section 3 for 
our dynamical results in the Conclusions. 
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Figure 5: The mao and fao in lattice units obtained from the fit of the dynamical scalar correlators 
with the PQChPT prediction (7) for various mo- The PQChPT prediction for two-flavor dynamical 
correlators incorporates oq and vrr/' exchange. The asterisks represent the fit for niQ oo, when the 
contribution from irrj' state vanishes and the fit reduces to the standard exponential fit (asterisks 
are the same as asterisks in Fig. 2). The different data points are slightly shifted from rUg = 
0.02, 0.03, 0.04 in horizontal direction for clarity. 



The extracted rriao and fao are almost insensitive to the presence of the tti]' state since 
the contribution of this state is at least an order of magnitude smaller than the dynamical 
lattice correlators in the fitted time-range for mo > 600 MeV. Our dynamical correlators 
are dominated by the oq exchange although the mass M^ + M^/ = M.„ + (M^ + ^mlY^'^ is 
comparable to the mass rriao in Table 2. This can be attributed to our spacial volume 16^ 
which is large enough to suppress the contribution e"^'^'^'^^'^''^'^^ /N'l of the irr]' exchange in 
Eq. (13). 



5.2 Analysis of the partially quenched correlators with rriyai 7^ i^sea 

Our analysis of the partially quenched correlators is based on the correlators with nisea = 0.02 
since a range of m^ai = 0.01 — 0.05 is available only in this case (see Table 1). 

The lattice correlators (1) and the PQChPT bubble contribution B{t) (7,8) are compared 
in Fig. 6 for m^ea = 0.02 and various myai- The PQChPT predictions in Figs. 6b and 6c 
show just the bubble contribution without the ag-exchange contribution in order to indicate 
the importance of the bubble contribution for various m^ai at fixed rrisea- Fig. 6 exhibits 
good qualitative agreement between the lattice correlators and the PQChPT predictions at 
various m^ai- 

• The lattice correlator and the bubble contribution are both negative and large for 
rriyai < rnsea due to a striking unphysical effect of partial quenching. The bubble 
contribution is large and negative for rriyai = 0.01 since it falls like — t-e"^^^^* at large 
t (9), where the corresponding pion mass Myv = 0.222(3) is small (Table 3). 
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Fi gure 6: The scalar correlator for nisea — 0.02 and various fTiyaV- the lattice data (a) and the 
bubble contribution B{t) as predicted by PQChPT for mo oo (b) and mo 800 MeV (c). The 
PQChPT prediction (7) in this figure represents just B(t) and does not contain the contribution 
from the ao-exchange. 



• The dynamical lattice correlator with mvai='^sea is positive. The bubble contribution 
describes the exchange of physical vrr/' and vanishes in the limit mo oo (12). 

• The lattice correlator at m^ai > ^sea is positive and receives a positive and rather small 
contribution from the bubble, which is less and less important for larger rriyai- The 
bubble contribution is positive and relatively small since it falls as a linear combination 
of -\-e~'^'^^st g^Yid +t ■ e~^^^^^* at large t (9), where the pseudoscalar masses Mys and 
Mvv are relatively large for rriyai > 0.03 and rrisea = 0.02. 

While we do not include it in the extraction of our final results, we have a limited set of 
partially quenched data for the rusea = 0.03 evolution, and we have checked that the lattice 
correlators change sign at rriyai = fTigea also in this case. This can be seen in Fig. lb. 

Let us have a closer look at the case of iriyai = 0.01 and mgea = 0.02 in Fig. 7, where the 
effect of partial quenching is most striking. The bubble contribution (dot-dashed) is in good 
quantitative agreement with the data for t > 8, where the ao-exchange fades exponentially. 
The exchange of the Cq scalar meson is dominant at smaller t. The fit of the data to the 
PQChPT prediction (7) gives m„o = 0.87(17) and /ao = 0.040(15) at mo ^ oo. The PQCHPT 
prediction with this choice of m^o and fao is given by the solid line in Fig. 7 and describes 
the lattice correlator well. Our one-loop analytical formula therefore correctly determines 
the sign and approximate size of the effects when the valence quark mass is lower than the 
sea quark mass. This gives us confidence in the veracity of applying this formula to the 
larger valence quark masses, where the loop effects are smaller. 
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t 

Figure 7: The data and PQChPT predictions for the scalar correlator at rusea = 0-02, m^al = 0-01 
and mo = oo. 

All scalar correlators for rUsea = 0.02 and various nivai arc fitted by the PQChPT pre- 
diction (7), and the resulting and fao are given in Fig. 8 and Table 4. Figures on the 
left represent the fit, which incorporates both the bubble and the ao-exchange contributions. 
We find that mao and /ao at m,^ai > 0.02 depend very slightly on the hairpin insertion mo 
in the range mo = [600 MeV, oo]. In the case of m^ai — 0.01, the central values of mao and 
/ao depend significantly on mo, but they are all consistent for mo in the range [600 MeV, oo] 
within the large error-bars^^. The result of the linear extrapolation from m^ai = 0.01 — 0.05 
to the chiral limit rriyai — > is practically independent of whether the rUyai = 0.01 data is 
taken into account due to the large error-bars at rriyai = 0.01 with current statistics. The 
linear extrapolation from rriyai — 0.01 — 0.05 gives 

mao = 0.90(11) , fao = 0.048(11) for mo oo 

mao = 0.89(9) , fao = 0.044(9) for mo = 800 MeV (15) 
mao = 0.88(9) , fao = 0.043(9) for mo = 600 MeV , 

which are consistent for mo = [600 MeV, oo] . So the chiral extrapolation m^ai — > at fixed 
'iTT'sea = 0.02 leads to the mass in the lattice units 

m^o = 0.89(11) , (16) 

where the error refiects the statistical error of the data and the variation of the bubble 
contribution for mo in the range mo = [600 MeV, oo]. 

^^Large error-bars on niao and fao at nival = 0.01 arise since they are oi)tained from the fit to the difference 
of the lattice correlator and the bubble contribution. The lattice correlator and the bubble contribution are 
negative and large for rrivai = 0.01, so their difference is small and has relatively large error-bar. 
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Figure 8: The niao and fao in lattice units obtained from the fit of the scalar correlators at 
iTT'sea = 0.02 with the PQChPT prediction (7). The left figures represent the fit results, when 
the bubble contribution is taken into account and tuq is varied. The right figures represent the 
conventional exponential fit e"'"""* + e~'""f(^^~*), which is obtained under the assumption that 
the bubble contribution vanishes; the correlator with m^ai = 0.01 and m^ea = 0.02 is negative and 
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horizontal direction for clarity. 
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The different data points are slightly shifted in 



The conventional exponential fit of the scalar correlators for rriyai > f^sea = 0.02 gives 
rriao and fao in Fig. 8 on the right. The bubble contribution in Eq. (7) is taken to be zero 
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in this case. The exponential fit obviously does not work for the correlator at rriyai = 0.01, 
where the intriguing partially quenched artifact has to be incorporated through the bubble 
contribution. However, it gives reasonable rriao and fao for m^ai ^ nT-sea'- the results from 
the exponential fit are consistent with the results from the fit to Eq. (7) since the bubble 
contribution is zero or relatively small for rriyai > ^sea- 

6 Non-singlet scalar meson mass 

In this section we collect our main results on the scalar meson mass. 

The chiral extrapolation of the two-flavor dyraamica/ data points myai—'>Tisea — '>TT'q 

gives (Eq. 2) 

ruao = 0.93 ± 0.20 or m^o^ = 1.58 ± 0.34 GeV , (17) 

where only the statistical error is given. The number in GeV is obtained using the preliminary 
result for the scale a"^ fsl.T GeV [14]. 

We extracted also the scalar meson masses from the partially quenched correlators with 
'f^vai^'f^sea- The chiral extrapolation m„ai — > at fixed rusea — 0.02 leads to 

m„o = 0.89 ±0.11 or mj^ = 1.51 ± 0.19 GeV (18) 

and we expect that the dependence on the sea quark mass is small. Here the error reflects 
the statistical error of the data and the variation of the bubble contribution for tjiq in the 
range itiq = [600 MeV, oo] (see Eq. 16 and Fig. 8). 

The chiral limits of Wao in the dynamical case and in the partially quenched case are 
consistent. Note that the error is smaller in the partially quenched case, where the application 
of the Partially Quenched ChPT was crucial. The mass of the simulated qq state is somewhat 
larger than the mass of ao(980) and it is closer to the mass of ao(1450). We note that 
our result is consistent with the fully quenched results of Refs. [8-10] within the present 
accuracy . 

Finally we comment on the scalar iso-doublet mesons su and sd, since their relation to the 
controversial resonance k is still an open question. We are not able to make a reliable estimate 
for the mass of the su and sd scalar mesons, since we did not simulate non-degenerate valence 
quarks. We get a rough estimate by extrapolating the mass obtained from the dynamical 
correlators to |(wis + w^M,d). The resulting mass ~ 0.92(9) or m^'*^ ~ 1.6 ±0.2 GeV seems 
higher than the mass of the reported experimental candidate ( ~ 800 MeV [5]) although the 
interpretation of this observation as k is controversial. 

^^The fully quenched results of [8 10] are presented in the Introduction and are consistent with (17,18) 
if the effect of quenching is incorporated at the leading order in the chiral expansion (one bubble). The 
quenched niao is somewhat lower if quenching effect is incorporated at the next-to- leading order [10]. 
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7 Conclusions 



Wc presented a lattice study of the lightest scalar qq state with non-singlet flavor structure 
(ao meson). Good chiral properties of the Domain Wall Fermions are important for the 
connected scalar correlator since this is the first step toward a controlled investigation of 
the scalar spectrum, in particular, the a particle, which is intimately related to the chiral 
symmetry breaking. Two degenerate sea-quarks were simulated with masses corresponding 
to ~ 500 — 700 MeV. The simulations were done at fixed lattice spacing and one size of 
the volume. 

The value of scalar mass = 0.93 ± 0.20 in the lattice units was extracted in the 
conventional way from the dynamical correlators (my^/ = m^e^) and the resulting error is 
rather large. The corresponding physical mass m^Q^ = 1.58 ± 0.34 GeV was obtained using 
the preliminary result for the scale a"^ ftil.7 GeV. 

We analyzed also the partially quenched correlators with rriyai 7^ i^sea- They exhibit 
striking effect of partial quenching since they are positive for m^ai ^ "^^sea and negative for 
T^vai < iTisea (Fig- !)• In Order to understand this effect of partial quenching, we derived the 
scalar correlator within the Partially Quenched ChPT. The leading unphysical contribution 
comes from the exchange of two pseudoscalar fields and has no unknown parameters. We have 
shown that this contribution is positive for m^ai > nT-sea, it is negative for rriyai < frisea and it is 
inversely proportional to spatial- volume at large time-separations. The physical contribution 
to the scalar correlator is due to the exchange of the scalar meson and has conventional 
exponential timc-dcpcndcncc. Wc find that the sum of these two contributions describes our 
partially quenched lattice correlators very well. We extract the mass = 0.89 ± 0.11 or 
m^o^ — 1.51±0.19 GeV from partially quenched correlators, which is consistent with the 
mass extracted from our dynamical correlators. 

Our current simulation of the qq state on the lattice seems to indicate that this state 
is somewhat heavier than the observed resonance ao(980), and it is closer to the observed 
resonance ao(1450); however, given the size of our errors this is not conclusive. We must 
also emphasize that we have only two dynamical flavors, our lattice volume is not very large 
and also our quark masses are quite heavy. Besides, continuum limit has not been taken as 
wc have data at only one lattice spacing, so the exploratory nature of our study needs to be 
kept in mind. 

We also note that the Nj = 2 theory is likely to have interesting differences from QCD 
{Nf = 2 + 1). Recall that the observed resonances ao(980) and ao(1450) decay to ?77r and 
KK. Bose statistics and isospin conservation restrict ao^~ 7^ tt -|- tt for Nj = 2, though 
aQ^~^ — > rj'^f=^ + TT would be possible if kinematics allows it. Additional intricacy could 
be also caused by the presence of a large 4-quark component in these channels. Thus the 
approach to the chiral limit may well exhibit a more involved dependence of the scalar mass 
on the quark mass than our data (sec Fig. 2) indicates with relatively heavy quarks. These 
issues will need to be addressed in future works with more computing resources. 
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Appendix A: Calculation of the bubble diagram 

In this appendix we derive the result B{p) (Eq. 8) for the bubble diagram in Fig. 3b. The 
coupling of the non-singlet scalar current and two pseudoscalar fields $ is given by Eq. (5) 

q2qi = 2//o(*')i2 ■ (19) 

The scalar correlator receives the following contribution from the exchange of the two pseu- 
doscalar fields shown in Fig. 3b 

B = {Q\q2qi qiq2\0)Fig.3b = 4/.^(0|$i„$„2 $26$w|0) , (20) 

where indices a and b are summed over all quarks and ghost-quarks of the partially quenched 
theory: a,b = S with i = 1, .., N^ai and S — 1, .., Nf. The non-zero Wick contractions 
relevant to the diagram on the Fig. 3b are 



+ ($ll|$22)(*12|*2l) + (*22|*ll)(*12|$2l) ) • (21) 



The propagators for the pseudoscalar fields follow from the Lagrangian (3). The propagator 
for the fiavor non-diagonal mesons (a ^ b) in Minkowski space is 

{^a.\^.a)=^J^, (22) 
y aa 

while for the diagonal mesons the propagator is [15] 

Sab^a 1 P^-Mls 1 



.P' - Ml, Nf (p2 - Mi) (p2 - Ml) 1 _ v^-Mh 



(23) 



The analytical expression (Eq. 8) for the bubble diagram in Fig. 3 is obtained by inserting 
the propagators (23) and (22) to the expression (21) and by performing the Wick rotation 
to the Euchdean space. 
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Appendix B: PQChPT correlator for a finite lattice 

The PQChPT prediction for the scalar correlator C^^'"^^^{t) (7) relevant for a finite lattice 
of the volume N^Nt is 

Nl/2~1 Nt /2-1 



ni,2,3=-iVi/2 n4=-NT/2 

1 1 + 

AT/ 



1 1 

+ 



with 

e = X:[2sin(i + [2sin(i (25) 

+ = X:[2sm(i If )]^ + [2sin(i{^ + ^})]^ ■ 



i=l 



The notation is given in section 4. This correlator is compared with the lattice correlators 
C^'^*(t) (6) in section 5. 
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rrivai 


configs. 


0.02 


0.01 - 0.05 


94 


0.03 


0.02 - 0.04 


94 


0.04 


0.04 


94 



Table 1: The summary of the scalar and pseudoscalar correlators analyzed in this work. All 
correlators are calculated atV = 16^ x 32, « 1.7 GeV, degenerate valence quarks, two degenerate 
sea quarks, point source and point sink. 





rriao 


faO 




t 


x' 


dof 


0.02 


0.92(9) 


0.044(9) 


4 


- 10 


0.1 


5 


0.03 


0.99(10) 


0.048(13) 


5 


- 10 


0.8 


4 


0.04 


0.94(5) 


0.042(6) 


5 


- 12 


0.2 


6 



Table 2: The niao and fao in lattice units obtained from the exponential fit to the dynamical 
correlators with m^ai = nT'sea- Time ranges t = tmin ~ tmaxt X and degrees of freedom {dof) in the 
fit are also shown. 
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x' 


dof 




0.01 


0.222(3) 


8-15 


0.4 


6 




0.02 


0.294(2) 


8-15 


0.06 
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0.02 


0.03 
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6 




0.04 


0.405(2) 


8-15 
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0.05 
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0.02 


0.304(2) 


8-15 


0.2 
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0.03 


0.03 


0.362(2) 


8-15 


0.05 
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0.04 


0.412(2) 


8-15 


0.02 


6 


0.01 


0.01 


0.108(2) 


8 - 15 


0.1 


6 



Table 3: The pion mass M^r in the lattice units, which is obtained from the fit to the pseudoscalar 
correlators at various mgea £^nd m^al- 
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Table 4: The m^o and fao in lattice units obtained from the fit of the correlators with the PQChPT 
prediction (7) at mgea = 0.02 and various ruyai- The fits denoted by "bubble and ao-exchange" 
take into account both terms in Eq. (7). The bubble contribution in (7) depends on the value 
of mo, which is taken to be niQ = 600 MeV, 800 MeV, oo. We also present the results of the 
conventional exponential fit e^™""* + e^™""''^^~*\ which is obtained under the assumption that the 
bubble contribution vanishes; the correlator with rrivai < fnsea is negative (Fig. 1) and can not be 
described by e"'"""* -|- e-'^''°^^'^~*\ 
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